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Although observations point to the neutrality and lack of currents on large scales in the universe, 
many mechanisms are known that can generate charges or currents during the early universe. We 
examine the question of survivability of relic charges and currents in a realistic model of the universe. 
We show that the dynamics of cosmological perturbations drive the universe to become electrically 
neutral and current-free to a high degree of accuracy on all scales, regardless of initial conditions. 
We find that charges are efficiently driven away in a time small compared to the Ifubble time for 
temperatures 100 GeV > 1 eV, while the same is true for currents at all temperatures T ^ 1 eV. 
The forced neutrality relaxes constraints on the generation of an electric charge in the early universe, 
while the forced erasure of currents disfavors many mechanisms for the early origins of large-scale 
magnetic fields. 

PACS numbers: 11.15.Ex, 95.30.Cq, 95.30.Qd, 98.62.Ai, 98.80.Cq 



The universe, on its largest scales, appears to be elec- 
trically neutral and current-free. Many scenarios exist, 
however, which can lead to the production of either a net 
electric charge or an electric current in the early universe. 
The presence of either charges or currents at early times 
has been discussed in the past as possible candidates for 
generating seed magnetic fields in the early universe, as 
well as causing other consequences for cosmology. In 
this letter, we expand upon our previous work [1], and 
demonstrate that although relic electric charges and cur- 
rents may reasonably exist at early times, they are driven 
away in a realistic, inhomogeneous universe. 

Many symmetries which appear to be good symmetries 
today may not have been so in the past. Such symme- 
tries are often related to conserved quantities, such as 
global symmetries (conserving baryon and lepton num- 
bers), or gauge symmetries (conserving color and electric 
charges). Although no direct experimental evidence ex- 
ists for the violation of these conservation laws [3| , the 
matter-antimatter asymmetry is a compelling indicator 
that perhaps baryon number conservation was violated 
at some point in thepast 0]. Indeed, so long as the 
Sakharov conditions [4| are met, baryon number viola- 
tion appears to be quite likely. Many mechanisms exist 
that could explain the generation of a baryon asymmetry 
(baryogenesis) , including the_decay of GUT-scale parti- 
cles , sphaleron processes [oj at the electroweak scale 
[tJ, sphaleron processing of a lepton asymmetry or 
through a coherent scalar field [91]. 

In a similar fashion to the global symmetry preserving 
baryon number, it may be possible to break the gauge 
symmetry preserving electric charge conservation. These 
scenarios can result in a process we refer to as electroge- 
nesis, since they admit the production of a net electric 
charge in the universe. Perhaps the simplest pathway 
for electrogenesis is a theory where the electromagnetic 
gauge symmetry U(l)em was temporarily broken in the 
past, only to be restored later at lower temperatures. 



Langacker and Pi [l^ demonstrated this possibility in the 
context of grand unification, while Orito and Yoshimura 
m and Nambu 12^ showed that electrogenesis is possi- 
ble in higher dimensional theories such as Kaluza-Klein. 
Another possibility that could lead to electrogenesis is 
that the U(l)cm symmetry is not exact. This arises in 
variable speed-of- light cosmologies [3], varying-a theo- 
ries [3], extensions of the standard model with massive 
photons ^Td\, brane-world models [l^. and models ad- 
mitting electron-positron oscillations [l7[. 

Additionally, electric currents can be generated at very 
early times by many means, such as during inflation [l8|, 
during the QCD phase transition 1^, 2ol| , the electroweak 
phase transition (2^, 21 1, or by relativistic decays [2^ . 



The earliest stages of structure formation also gener- 
ate small electric currents |H, but these are too small 
to be of any significance until much lower temperatures 
(T <C 1 MeV) . Prior treatments of a universe with some 
sort of electromagnetic asymmetry have assumed that 
charges and/or currents within a comoving volume re- 



main constant |lll . |22| . |23| . |2J]. However, charge con- 



servation still allows the local charge to change, provided 
there is a flow of current, and therefore there can be local 
changes in the charge-per-baryon (A) at different epochs. 
Observational constraints on A at different epochs are 
as follows: the anisotropics of cosmic ray s place a con- 



straint that at present, |A| < 10 e [llj; the isotropy 



of the microwave background constrains | A| < 10 ^ e at 
z ~ 1089 25]; and big bang nucleosynthesis requires that 
]A] < 10-32 e at z ~ 4 X 10^ [i^. 

In order for a local cosmological charge or current to 
have any major signiflcance, they must be correlated 
with the density inhomogeneities which will grow into 
the collapsed structure observed today. It is expected 
that charge over- and under-densities should possess the 
same types of inhomogeneities as baryons ,22], and their 
evolution should therefore be calculable in the same fash- 
ion. We have recently shown that gravitational forces, in 
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combination with Coulomb forces and Thomson scatter- 
ing, all impact the evolution of a charge asymmetry [l[ . 
The remainder of this letter examines the evolution of a 
relic electric charge or current, and studies the associated 
cosmological ramifications. 

A realistic cosmological model of the universe will nec- 
essarily include inhomogeneities on both subhorizon and 
superhorizon scales, as mandated by inflation. An excel- 
lent treatment of the linear evolution of these perturba- 
tions, including the photon, neutrino, baryonic, and dark 
matter components, is given in Ma and Bertschinger 27 1 . 
While it is possible that a cosmological charge asymmetry 
could manifest itself in exotic forms, these possibilities 
are unsupported by experiment. This includes searches 
for fractionally charged particles 24, 2^, or a difference 



between the strength of the proton and electron charges 
[2^ . We therefore consider that if a cosmic charge asym- 
metry exists, it is due to a difference between the number 
densities of_protons and electrons, following our previous 
treatment [l]. 

We choose to work in the conformal Newtonian gauge, 
defined by the metric 

ds^ = a2(T)[-(l + 2i^)dT^ + (1 - 2(j))dx'dx^], (1) 

where (p and ip describe the scalar-mode inhomogeneities. 
Including Coulomb interactions and Thomson/Compton 
scattering with photons, the evolution equations for the 
proton and electron fluids become 



S, = -(. 



30, 



0, + cies, + fc^V 



(2) 
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where Si is deflned as the the departure of the local spa- 
tial density of species i from the overall spatial average 
(i.e.. Si = pi(x)lpi — l), and 0i is defined by 0i = iy{vi)j, 
where Vi is the peculiar velocity of the field 5i. In equa- 
tion ([2]), i is p for protons and e for electrons, rip and 
rie are the local number density of protons and electrons, 
qp and (/e are the proton (-l-e) and electron (— e) charges, 
and Ti is the (conformal time) rate of momentum transfer 
due to photon scattering with charged particles. 
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The Thomson cross section (ctt) for electrons is replaced 
with the Klein-Nishina form for temperatures T ^ me 
[3^. By taking the difference between the density and 
velocity fields for electrons and protons, a set of evolution 
equations governing the evolution of a net charge and/or 
current within a given volume is obtained. The equations 



for Sq and 9q, where dg 



(4) 



9q + cik'Sq 



%+0q)+LU^Sq, 



to linear order, where = {mp5p -f rrieSe) / {mp + rrie) 
is the baryonic mass density perturbation and uj"^ — 
A-nriee^a? /me is the (conformal time) plasma frequency. 
Note that other types of scattering, such as Rutherford 
scattering, do not need to be included in equation (j4]) 
at this epoch, as the momentum transfer arising from 
Thomson/Compton scattering is by far the dominant 
scattering term in the early universe. 

In our previous paper jlj , we investigated the effects of 
the source term in equation proportional to 0^ — 0b, 
in the absence of any initial charge or current asymme- 
try. We found that a local charge asymmetry, of order 
6q ^ 10^^'* on Mpc scales, is generated near the time of 
recombination. In this paper, we study the fate of an 
initial charge (Sq) or current (6q) asymmetry. We note 
that since equations ([4]) are linear, these two cases evolve 
independently of one another. As the rates Te and lo 
are typically much greater than the expansion rate d/a 
before decoupling, and much greater than the term c^fc^ 
on all but the smallest scales, the latter terms can be 
neglected, and equations (jH) can be simplified into the 
single equation 



TeS. 



=0, 



(5) 



where we remind the reader that 5q = —Oq. An initial 
asymmetry evolves as a damped harmonic oscillator, with 
slowly varying coefficients. If the coefficients were con- 
stant, we would have the usual solution, with two modes. 



5q = Cie*'^ -I- 026^* 



(6) 



where ci and C2 are undetermined constants, and si,2 are 
the two roots of — VeS + cj^ = 0, 



Sl,2 = -Te/2±^Tl/^-L0\ 



(7) 



For Te/2 > u!, the modes are a fast and a slow decay, 
while for Te/2 < uj, there is an oscillation with a damped 
envelope. 

The rates Te and uj in fact vary with expansion, as 
Te ex. and uj oc a}/"^, and so change slowly on the 
expansion timescale. As a result, the arguments of the 
exponentials in equation Q change from s t to J s dr, as 
can be seen straightforwardly by changing variables to 
the logarithm of the asymmetry, X — — ln(5g/(5o)j where 
(5o is the initial asymmetry, and its derivative, the damp- 
ing rate s — 9q/5q. From equations these quantities 
satisfy 



X = s. 



s = s 

a 



+ cik^ + CJ^ - TeS + , 



(8) 



where the rates uj and Te are again much greater than the 
expansion rate H = d/a and the term c^fc^ on cosmolog- 
ical scales. We expect s is also small, since the rates Te 
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FIG. 1: The logarithm of the damping factor X = — X-aiiqj&d) 
(solid curve) and the normalized damping rate per Hubble 
time sjH (dashed curves), as a function of redshift. The lower 
dashed curve represents the slow decay mode, which describes 
the evolution of an initial charge, while the upper dashed 
curve illustrates the fast mode, describing the damping of 
an initial current. After critical damping (z ~ lO'"*), the two 
dashed curves evolve together. For damping rates s/H > Iff , 
any initial charge or current will be completely driven away in 
less than one Hubble time, which happens for excess charges 
for values of z between 10'^ and 10^"^ and for currents for 

z > lo^ 



and to vary slowly. In this approximation, the damping 



rate s is a root of s — r„s 



0, as above, and 



X 



i dr = 



s f da s 

da - ^ / — 77- 

a J a H 



(9) 



The quantity s/H represents the rate of logarithmic 
damping per Hubble time, or damping per log expansion 
factor. 

Figure [T] demonstrates dramatically that the damping 
is so powerful as to wipe out a charge or current asym- 
metry for a large window around the epoch of critical 
damping. The figure shows the logarithm of the damp- 
ing factor for both a charge asymmetry 5q oc e~'^ and a 
current J (x 9 oc se~'^ . At early times, an initial charge 
asymmetry decays according to the slow-decay mode so- 
lution of equations ([HI), while an initial current follows 
the fast-decay mode. For values of X > 10^, an initial 
charge asymmetry of order 1 is reduced to less than one 
excess proton or electron per comoving horizon. 

The results in Figure [U indicate that any excess charge 
created at temperatures above 100 GeV will be entirely 
wiped out at roughly the electroweak scale, and any elec- 
trogenesis that concludes between Tew and Tcmb will 
also be driven away in a time much less than the Hubble 
time. Furthermore, any relic current created at temper- 
atures T 3> Tcmb will be also driven away in a time 
much less than the Hubble time. We therefore conclude 
that the evolution of charged particles in the universe, in- 



cluding Thomson/Compton scattering and the Coulomb 
interaction, force the universe to an electrically neutral, 
current-free state, independent of initial conditions, infla- 
tionary remnants, electrogenesis, and phase transitions 
that may leave charges or currents behind. 

The results presented in this paper may have signif- 
icant implications for scenarios admitting electrogenesis 
in the early universe. The observational limits on the 
charge-per-baryon (A) of the universe H, 25, [il] have 
been used to constrain models admitting electrogenesis in 
the past. Our results demonstrate that later-time con- 
straints on A do not constrain initial values of A, as the 
universe drives away any local charge asymmetries expo- 
nentially quickly. Therefore, we also conclude that the 
U(l)em gauge symmetry may have been broken in the 
past or may be only approximately exact today, and the 
observed values of A would not be discernably different 
from zero. Hence, electrogenesis is not constrained at 
all by observational limits on A. Nevertheless, the phe- 
nomenological consequences of electrogenesis may leave 
signatures which can be searches for experimentally. Just 
as GUT-scale baryogenesis has phenomenological conse- 
quences such as proton decay, electrogenesis should also 
have consequences for particle physics. Each specific 
model that admits electrogenesis will have its own as- 
sociated phenomenology. Some consequences of models 
admitting charge nonconservation have been worked out 
[30| . with various distinct experimental signatures hav- 
ing been (unsuccessfully) searched for |3lj] . The physical 
possibilities if either Lorentz invariance or the U(l)cm 
symmetry are broken is quite rich and varied, and de- 
serve further examination. 

Also of cosmological importance are the implications 
for magnetic fields in the early universe. If large-scale 
currents are truly driven away, as in Figure [U then not 
only are magnetic fields generated by a charge asymmetry 
[23] ruled out, but all magnetogenesis scenarios which 
rely on a relic current persisting from early times are 
disfavored as well. This requires careful examination, as 
the scenarios for magnetic field generation from the QCD 
[3,12^, electroweak [2^, 21|, or inflationary 18 1 phase 
transitions correctly discuss the growth of instabilities 
in the hot plasma of the early universe. However, they 
do not examine whether the currents which are initially 
generated by the phase transition persist long enough 
to create such magnetic fields. Our results in Figure [T] 
illustrate the damping rate for currents per Hubble time, 
and we find the damping rate to be much quicker than the 
growth rate of the instabilities given in the literature. If 
the currents do not persist long enough to create small- 
scale fields, they will have no opportunity to then be 
transferred to larger scales [S^l, creating seeds for the 
large-scale magnetic fields observed today. Since these 
currents do not have a source driving them after their 
creation, the expansion dynamics of the universe ought 
to damp them away very rapidly. Unless there is some 
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additional effect which prevents the currents that would 
give rise to magnetic fields from being driven away, these 
magnetogenesis scenarios are highly disfavored. 

Since our results indicate that initial conditions on Sq 
and 0q are unimportant in the absence of sources, the 
most significant contribution to a late-time Sq or 9q must 
come from the source contribution in equations ^ . After 
any initial charge is wiped out, Sq and 9q evolve according 
to the dynamics in equations ([4]), where the source term, 
re(^7 — db), determines Sq. Results from our previous 
work [l[ indicate that the typical charge-per-baryon on 
the scale of the present horizon is | A| ~ 10~^^ e, a result 
that ought to be valid since the present horizon is still in 
the linear regime of structure formation. Strictly, these 
calculations apply to asymmetries within the horizon and 
at finite Fourier wavenumber k. The asymmetry must be 
within the horizon for its dynamics to evolve according 
to our analysis due to causality; thus our results are defi- 
nitely valid for all wavenumbers k from the moment they 
first enter the horizon. We also consider a naive (and not 
rigorous) extrapolation to superhorizon ks. The depen- 
dence of our results on k is unimportant and vanishes as 
A; — > 0. Thus, it appears that even a global charge asym- 
metry may be driven to zero. At first glance, this may 
seem to violate the principle of charge conservation, but 
we recall that currents carry away excess charge while 
still satisfying the continuity equation. An intuitive pic- 
ture is to consider an expanding universe containing ho- 
mogeneous proton and electron fluids free to expand at 
different rates. In the Newtonian limit, with spherical 
symmetry, an excess of positive charge increases the ex- 
pansion rate for protons and decreases the expansion rate 
for electrons, and the total charge in a given volume can 
decrease even though charge is locally conserved. This 
construct does not work in detail, as different expansion 
rates require a radial flow of protons relative to electrons, 
and such a universe could be isotropic about only one 
center. The flow may also require relativistic velocities. 
However, the mechanism can apply on arbitrarily large 
finite scales, and to reduce a charge asymmetry of order 
the mass density inhomogeneity (10~^) over the scale of 
the horizon requires a relative velocity of 0(km/s), or a 
bulk fluid displacement of < 1 Mpc. 

Finally, although our calculations have accounted for 
high-energy effects in electron-photon scattering [s^l , we 
have not done this for protons. We note that protons 
were dissociated into a quark-gluon plasma at tempera- 
tures above the Aqcd scale, resulting in different scat- 
tering rates and momentum transfer at these energies. 
Furthermore, we have not taken into account possible 
shielding effects on cosmological scales, which, if present, 
could cause the effective Coulomb force to fall off faster 
than l/r^ for sufhciently large distances. Inclusion of 
these effects may be needed to extract the exact behav- 
ioral details of Sq and dq at high energies. 

The overall conclusion which can be drawn from this 



work is that any relic charge or current created on sub- 
horizon scales in the early universe (at temperatures 
T ^ 1 eV) is driven away by cosmological dynamics. 
This indicates that a net charge or current of any magni- 
tude could be generated at the QCD phase transition, the 
electroweak scale, the time of Higgs symmetry breaking, 
the time of supersymmetry breaking, the end of inflation, 
or the grand unification scale, and the universe would not 
be discernibly different from a universe that was electri- 
cally neutral and current-free at all times. This disfa- 
vors many scenarios for the generation of seed magnetic 
fields, which rely on the survival and sustainability of 
currents at early times. The possibility that the universe 
underwent electrogenesis at some point is intriguing and 
rich with possibility, and present observations of electric 
charge neutrality do not preclude the possibility that the 
universe (or a portion of it) was charged in the past. Al- 
though the phenomenological consequences of electrogen- 
esis have been elusive thus far, evidence for its existence 
would profoundly alter our cosmological picture of the 
universe. 
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